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ABSTRACT 

Net tidal torque by the secondary on a misaligned accretion disk, like the net tidal 
torque by the Moon and the Sun on the equatorial bulge of the spinning and tilted Earth, 
is suggested by others to be a source to retrograde precession in non-magnetic, accreting 
Cataclysmic Variable (CV) Dwarf Novae systems that show negative superhumps in 
their light curves. We investigate this idea in this work. 

We generate a generic theoretical expression for retrograde precession in spinning 
disks that are misaligned with the orbital plane. Our generic theoretical expression 
matches that which describes the retrograde precession of Earths' equinoxes. By mak- 
ing appropriate assumptions, we reduce our generic theoretical expression to those gen- 
erated by others, or to those used by others, to describe retrograde precession in pro- 
tostellar, protoplanetary, X-ray binary, non-magnetic CV DN, quasar and black hole 
systems. 

We find that spinning, tilted CV DN systems cannot be described by a precessing 
ring or by a precessing rigid disk. We find that differential rotation and effects on 
the disk by the accretion stream must be addressed. Our analysis indicates that the 
best description of a retrogradely precessing spinning, tilted, CV DN accretion disk is 
a differentially rotating, tilted disk with an attached rotating, tilted ring located near 
the innermost disk annuli. In agreement with observations and numerical simulations 
by others, we find that our numerically simulated CV DN accretion disks retrogradely 
precess as a unit. Our final, reduced expression for retrograde precession agrees well 
with our numerical simulation results and with selective observational systems that 
seem to have main sequence secondaries. 

Our results suggest that a major source to retrograde precession is tidal torques 
like that by the Moon and the Sun on the Earth. In addition, these tidal torques 
should be common to a variety of systems where one member is spinning and tilted, 
regardless if accretion disks are present or not. Our results suggest that the accretion 
disk's geometric shape directly affects the disk's precession rate. 



Subject headings: accretion, accretion disks; methods: analytical; binaries: general; 
binaries: close; stars: dwarf novae; novae, cataclysmic variables 
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1. Introduction 

Variability in non-magnetic, Cataclysmic Variable (CV) close binaries has been largely at- 
tributed to instabilities. Outbursts in recurrent and classical novae can be caused by thermonu- 
clear runaway on the surface of a white dwarf or by a thermal and viscous instability in the disk. 
The latter outburst involves cycling from low-to-high thermal and viscous states (see e.g., Meyer 
& Meyer-Hofmeister 1981 and references within), and thus cycling from low-to-high mass transfer 
rates (see e.g., Lasota 2001), that results in a release of accretion energy as mass is drained onto 
the white dwarf. This type of outburst is similar to that experienced by CV Dwarf Novae (DN), 
but larger in scale, and by transient X-ray binaries that have neutron stars or black holes as the 
primary star in the binary. Another type of instability, a tidal instability, results in CV DN SU 
UMa and Nova-Like (NL) outbursts, of which the larger amplitude and longer duration superout- 
bursts involve both a thermal- viscous instability and a tidal instability. Both CV DN SU UMa and 
NL outbursts and superoutbursts show characterisitc hump-shaped modulations, known as positive 
superhumps, in their light curves. Positive superhumps have periods that are a few percent longer 
than the orbital period. 

In addition to positive superhumps, some non-magnetic CV DN systems (e.g., TT Ari) show 
negative superhumps in their light curves. These negative superhumps have a period that is a few 
percent shorter than the orbital period. Negative superhumps are suggested to be a consequence of 
a partially tilted disk (Patterson et al. 1993) or are related to a warped accretion disk (Petterson 
1977; Murray & Armitage 1998; Terquem & Papaloizou 2000; Murray et al. 2002; Foulkes, Haswell, 
& Murray 2006). Our previous numerical simulations (Wood, Montgomery, & Simpson 2000; 
Montgomery 2004, Montgomery 2009) show that negative superhumps can be produced by a disk 
fully tilted out of the orbital plane. In Montgomery (2009), we suggest that the source to the 
negative superhump is extra light from inner annuli of a tilted accretion disk, and this extra light 
waxes and wanes with the amount of gas stream overflow and particle migration to inner annuli as 
the secondary orbits. 

Accretion disks have been postulated to warp or tilt via a potpourri of sources. For example, 
disk tilt in X-ray binaries can be from gas streaming at an upward angle from the inner Lagrange 
point for one half of the orbit and at a downward angle for the second half of the orbit (Boynton et 
al. 1980). In some CVs, a disk tilt can be held constant by a gas stream that is fed via the magnetic 
field of the secondary (Barrett, O'Donoghue, & Warner 1988). Also for CV systems, a disk tilt 
instability can result from a coupling of an eccentric instability to Lindblad resonances (Lubow 
1992). A vertical resonant oscillation of the disk midplane can be caused by tidal interactions 
between a massive secondary and a coplanar primary (Lubow & Pringle 1993). A warping instability 
can be caused by irradiation from the primary (Pringle 1996, 1997). A warping can be caused by 
direct tidal forces from a secondary orbiting on an inclined orbit (e.g., Papaloizou & Terquem 1995, 
Larwood et al. 1996, Larwood 1997, Larwood & Papaloizou 1997). A disk tilt can be induced if 
the secondary has a body axis that is tilted relative to the orbital plane (Roberts 1974). A disk 
warp can be caused by misalignments of the spin axis of a compact and/or magnetized primary 
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and the disk axis (see e.g., Kumar 1986, 1989). To date, no consensus has been found regarding 
the source to a disk tilt or warp. We do acknowledge that Murray &; Armitage (1998) find that CV 
DN accretion disks do not tilt significantly out of the orbital plane by instabilities. 

Regardless of how an accretion disk tilts or warps, non-coplaner accretion disks should process 
in the retrograde direction, that is, the direction opposite to spin and orbital motion. Precession 
is the motion of a body axis around a fixed axis in space. Retrograde precessional motion traces 
out a cone around the cone's symmetric axis, a fixed axis in space. Two sources are known to 
generate retrograde precession. One is torques acting on the body. A second kind is fast precession 
which applies to bodies that are free. An example of the latter is a coin thrown into the air 
in such a way that the coin spins fast, but the normal to the coin wobbles. Of the two, many 
typically cite torques as the source to retrograde precession. For example, retrograde precessional 
activity has been suggested for X-ray binaries (e.g., Katz 1973, Larwood 1998, Wijers & Pringle 
1999), for common systems that produce jets (e.g., Livio 1999), for protostars (e.g., Papaloizou & 
Terquem 1995, Larwood et. al. 1996, Larwood 1997), for quasars and supermassive black holes 
(e.g., Romero et al. 2000, Caproni & Abraham 2002), for CVs (Barrett et al. 1988, Patterson 
et al. 1993, Harvey et al. 1995), to name a few. Katz (1973) suggests that the permanently 
tilted edge of a disk precesses due to the tidal field of the secondary. Roberts (1974) suggests 
that the secondary precesses due to the tidal field of the primary and the tilted disk retrogradely 
precesses as a consequence. Barrett et al. (1988) and Kondo et al. (1983) suggest that particles 
may be orbiting a tilted disk and as a result, the disk retrograde precesses. Patterson et al. (1993) 
notice a near 2:1 ratio of the negative-to-positive precessional periods of several CV systems that 
exhibit both modulations in their light curves, and they note that this ratio is curiously similar to 
that of the Earth which has a retrograde nodding precessional period P r =18.61 yr and a prograde 
precessional period P p =8.85 yr (see e.g., Stacy 1977). Patterson et al. (1993) suggest that the 
source to retrograde precession in these CV DN systems may be the same as that which causes the 
Earth to retrogradely precess. Petterson (1977) suggests that the Earth's retrograde precession is 
caused by outside forces, namely tidal torques on the Earth's bulge by the Moon and Sun with the 
Moon having the most effect due to its closer proximity to Earth. Hence, Patterson et al. (1993) 
suggests that a tidal torque by the secondary on a tilted accretion disk could be the source to 
retrograde precession in CV DN non-magnetic, accreting disks. 

Others have found that only portions of their disk retrogradely precesses. For semidetached 
binaries, Bisikalo et al. (2004) show that a one-armed spiral shock wave can be generated in inner 
annuli of cooler disks (i.e., T ~ 10 4 i"Q. They find that their inner annuli spiral wave precesses 
retrogradely, but the disk does not. Based on their numerical simulations, they show that the 
wave is generated by sheared elliptical orbits. They suggest that the sheared elliptical orbits are 
due to retrograde precession within the disk, and the retrograde precession is due to an influence 
by the mass-losing star by the processes discussed in Kumar (1986) and Warner (1995, 2003). 
However, this suggestion cannot be correct because Kumar (1986) finds that retrograde precession 
is due to the tidal influence of the secondary on a tilted disk's outer annulus, not inner annulus, 
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and because the simulations by Bisikalo et al. (2004) do not have warped or tilted disks (private 
communication) . 

In this work, we test the idea of Patterson et al. (1993) on the tilted accretion. We assume 
the disk is fully tilted out of the orbital plane because our previous works (Wood, Montgomery, 
& Simpson 2000; Montgomery 2004, Montgomery 2009) are successful in showing that fully tilted 
disks retrogradely precess. In addition, our works also generate negative superhump modulations 
in our artificial light curves that have the correct shape and negative superhump periods of the 
right percentage as compared with observations. In this work we do not explain how accretion 
disks tilt but instead check to see if the source to retrograde precession is like that which causes the 
spinning, tilted Earth to retrogradely precess. In §2 of this paper, we develop a generic theoretical 
expression to describe retrograde precession using the Earth-Moon-Sun system as a model. In 
§3, we reduce the generic theoretical expression by considering the geometry of the system and 
by making appropriate assumptions, thereby obtaining expressions generated or used by others to 
describe retrograde precession in protoplanetary, protostellar, X-ray, black hole, and non-magnetic 
accreting CV binary systems. In §4, we summarize our numerical simulations. In §5, we list 
observational data. Is §6, we compare retrograde precessional values obtained from our theoretical 
expressions with those obtained from our numerical simulations and with those obtained from 
observational data. In §7, we provide a discussion and in §8, we summarize, conclude, and list 
future work. 



2. Theoretical Expressions for Retrograde Precessing 

2.1. Binary Systems 

We consider two systems, one being the Earth-Moon-Sun system and one being a non-magnetic, 
accreting CV DN system where the accretion disk is misaligned with the orbital plane. In the 
Earth-Moon-Sun system, the spinning Earth is known to have an oblate spheroidal shape. The 
resultant equatorial bulge is not symmetric about Earth's rotation axis. As such, a net torque on 
this equatorial bulge results in the rotational axis of the Earth precessing in the opposite sense to 
the Earth's spin and Earth's orbital directions. 

To simplify the many-body problem of the Earth- Moon-Sun system and to generate like coor- 
dinates among the two systems, we only consider net torque effects by the Sun on the Earth. Figure 
1, top panel, exaggerates the Earth and Sun size and distance scales to show details. In this panel, 
the Sun is the small sphere on the y axis located at a distance d from the origin. The Earth is the 
large oblate spheroid that is centered at the origin in the (x, y, z) coordinate system of the ecliptic. 
The z axis defines the normal to the ecliptic and orbital motion Co is prograde. The (x',y',z') 
coordinate system (y' not shown) is the equatorial system of the Earth. Earth rotates around the 
z' axis with angular velocity e. As shown, the angle from the z to the z' axes is the obliquity 6. 
Note that x and x' share the same axes and thus they coincide with the line of equinoxes at this 
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instant in time in the Earth-Sun system (or with the line of nodes in the non-magnetic, accreting 
CV DN system). Also note that z', z, and y are in the same plane and thus Eulerean rotation is 
about the x, x' axis. Because the gravitational attraction forces differ from a point on Earth nearer 
to the Sun, denoted by the force vector F n , than a point on Earth further from the Sun, denoted 
by the force vector Ff , torques around the origin result. 

Figure 1, bottom panel, represents the CV DN system and is similar to Figure la except that 
the oblate spheroidal Earth is replaced with a massive primary star Mi that has angular velocity 
e and that is surrounded by a disk of mass Md where Md <C M\. Also, the Sun is replaced with 
a secondary star of mass Mi- In this non-magnetic CV DN system, the massive primary cannot 
be seen as it is smaller in size compared to the disk size. We assume that the separation d of the 
system is much larger than the thickness of the disk, and we assume that the disk has negligible 
mass compared to that of the primary. Like the top panel, this bottom panel is not to scale. 

Ultimately we seek the retrograde precession of the line of nodes in the spinning, tilted CV DN 
system. To find this precession, we use the precession of the equinoxes in the Earth-Sun system as 
a model. 



2.2. Angular Momentum and Moment of Inertia Tensor 

For continuous bodies, its angular momentum vector L is related to its angular velocity vector 
fi in any coordinate system through a linear transformation matrix known as the moment of inertia 
tensor I, 

l = m (i) 

Using a right-hand coordinate system whose general coordinates are (x',y',z'), we can identify the 
diagonal elements of the matrix as the moment of inertia coefficients and the off-diagonal elements 
as the products of inertia. If we assume, for example, that fl is about the x' axis and therefore L 
is along the x' axis, then 

L x ' = I x ' x 'Q x r + I x t y iVLyi + I x r z /Q z i. (2) 



If the density p of a solid is known, then the moment of inertia can be found from / = 
J p r /r' 2 dV' where r'j_ is the perpendicular radial distance from any rotational axis, dV' is small 
volume element, and density is a function of radius. In cartesian coordinates, the moment of 
inertia tensor is 



j P(x',y',z') 



y !2 + 

—x'y 1 
-x'z' 



—x'y' 
z' 2 + x' 2 
-y'z' 



—x'z' 
-y'z' 
x' 2 + y' 2 



dV 



(3) 



where dV'=dx'dy'dz'. If cylindrical coordinates are desired, then we can integrate from (—h c /2) 
to (h c /2) along the z' axis where h c is the height of the cylinder, from — to along the y' 

axis where is the radius of the cylinder, and from y' 2 to y^r J — y' 2 along the x' axis. 
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However, the solution to this integral is most easily found using curvilinear coordinates q± = r', 
Q2 = 4>, Q3 = z' with scale factors h\ = 1, hi = r', and h% = 1 for a circular cylinder. In this 
solution, dV = r' dr' d<pdz\ x' = r' cos <f), y' = r'sin</>, and z' = z' . The limits of integration arc 
— h c /2 — ► h c /2 on dz' , — > 2tt on g?0, and — ► on dr'. 



2.3. Potential 

In the Earth-Sun system, the potential we seek is that of the Earth and thus we center this 
attracting body at the origin of our coordinate system. As shown in Figure 2, the attracting 
body has mass m; an elemental mass dm is located at a distance r' from the center of mass of 
the attracting body and at an angle a to the line of centers connecting the center of mass of the 
attracting body to the center of mass of a far-away point mass object M. We assume in the Earth- 
Sun system that the Sun is located at a distance d where d» r'. The potential V at the far away 
point mass object is 



GMm GM 

V= -^3- (i*V + ly'y' + Iz'z> ~ 3/) (4) 



where we have assumed that Earth has uniform density p and rotation is possible about the line of 
centers. In this equation, G is the universal gravitational constant. 

This potential is that of MacCullagh, a result that is applicable for slightly non-spherical 
celestial bodies and for d » r'. This potential is valid for the oblate spheroidal Earth in the 
Earth-Moon-Sun system. We assume that this potential is also valid for the spinning, tilted CV 
DN system even though the primary is surrounded by a disk (and hence the geometry is by far 
not spherical) and the radius of the white dwarf primary plus accretion disk is not significantly less 
than the distance to the point mass secondary star. As the first term in the potential dominates 
over the other terms (of which only the second term is shown as all others are assumed negligible) 
and as <C Mi, the potential is most affected by M and m = Mi + M^ ~ Mi. As Mi is assumed 
spherical and of uniform density and as the secondary is significantly further away than compared 
with the radius of the primary, we accept this potential as valid for the spinning, tilted CV DN 
system. We note that the second term in Equation (4) is not negligible as the disk's geometric 
shape affects the potential via the moment of inertia. We also note that as the radius of the disk 
is dependent on the mass ratio of the system, with smaller mass ratios yielding larger disk radii, 
accretion disks in small mass ratio systems may be subject to additional perturbations. These 
additional correction terms are not considered in this work. 
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2.4. Net Torque 

In Figure 2, the Earth in the Earth-Sun system and the primary in the CV DN system are 
the attracting mass bodies that are centered at the origin. The Sun and the secondary star are 
considered to be the far away point mass objects. As the potential is evaluated at the far away 
point mass object, the force, which is proportional to the gradient of the potential, is on the far 
away point mass object. By Newton's Third Law, the far away point mass object exerts an equal 
and opposite force on the attracting body. As the forces are not co-linear, a torque about the origin 
arises. Figure 1 is exaggerated to show non-colinear forces. 

Because the spinning Earth has rotational flattening and because the Earth's spin axis is tilted, 
the force pulling on the near-side of Earth (i.e., the side closer to the Sun) is greater than the force 
pulling on the far-side of the Earth, and a net non-zero torque results about the x axis (see Figure 
1, top panel). In the spinning, tilted CV DN system, a torque about the origin arises for the 
same reasons - because the disk is tilted and spinning and because this combination results in 
non-colinear near and far side forces on the disk (see Figure 1, bottom panel). As shown in Figure 
1, the body is rotating counter-clockwise about the z' body axis, and the body moment of inertia 
coefficients are along the body axes (x',y',z'). 

By definition, the net torque around the origin of the body is T = r' x F where any force 
F = — V(V). As r' x r' = 0, then the I in Equation (4) is the only term that is not just a function 
of radius but also a function of (x',y',z'). Therefore, we reduce Equation (4) to 

3GM ,2 ,2 T i2\ fa\ 

= -^-{Ix'x'X +Iy>y>y +h'z'Z ) (6) 

where we have evaluated the potential at a distance r' = d from the origin. The net torque in the 
(x',y',z') coordinate system is 

r = r x , + r y , + r z , (7) 

—3GM \/ / /\ / 

= d5 (Iz>z> - ly'y'KV Z )X - 

3GM , , , 

— ^5— {I x 'x> - Iz>z>){x z )y - 

(-W ~ Ix'x'){x'y')z'. (8) 

By employing transformation equations (see e.g., Goldstein 1980), we can rotate the net torque 
from the (x',y',z') coordinate system to that of the (x, y,z) coordinate system (see Figure 1). In 
this transformation, we assume that the obliquity angle 9 remains constant, the orbit is circular, 
and the mean orbital motion is u> as shown in Figure 1. The resultant net torque in the (x, y,z) 
coordinate system becomes r = T x + T y + T z with component magnitudes 

3GM 

T x = ^ 3 {h'z' ~ Iy' y ')[sm9cos6sm 2 (ibt + ip)} (9) 
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3GM 

T y = ^ 3 (I X 'x' — Iz'z') [cos(cjt + if)) sm(ut + ip) sin 0] (10) 



r z = 



-3GM 



{ly'y' ~ I x ' x ')[cOs(d)t + if)) Sm(d)t + if)) COS 9]. (11) 



where ip is some phase angle, and t is time, and to is an angle that increases in the orbital direction 
from perihelion. In these equations, we have neglected any variations over one orbit because the 
precessions we seek are over long periods of time, much longer than one orbital period. 

To simplify the equations further, we assume symmetry about the z' rotation axis, and we 
assume a circular disk. Therefore, I x i x i = I y i y i and T z is zero. We cannot simplify T y any further 
but we we can substitute the identity sin(2#) = 2 sin 9 cos 9 into T x . Because the magnitude of this 
torque is not constant, we need to take averages. For example, the magnitude is a minimum at 
either of the equinoxes and a maximum at either of the solstices in the Earth-Sun system. Likewise, 
the magnitude is a minimum when the primary, secondary, and the disk's line of nodes align and 
a maximum when the primary, secondary, and disk's line of antinodes align. Therefore, we can 
find the average of T x and T y over a time period of one orbit using an integral average of some 
continuous function f(t) over an interval r, < /(<) >= \ Jq f(t)dt. If we assume in J^ x that only 
the sin 2 (cut + ip) is not constant and we assume the phase constant is zero, then we find that 
< sin 2 (a>i + ip) >= 1/2. For the same time interval, we find that T y = 0. Of Equations (9)-(ll), 
only T x is left, 

3GM 

Tx = ~4d3~ (4 '*' " Ix ' x,) Sin(20) ' (12) 
We can use unit vector notation and employ the cross product definition, 

z'xz = (z'z sin 0)x' 

= (sin0)x', (13) 

to obtain the component of the net torque along x, 

~3GM 



2d 3 



(I z > z > - I x > x >) cos(0) 



(14) 



For the Earth-Sun system, M = M sun and is not in solar masses. For the Earth-Sun system, the 
component of the average torque along the line of equinoxes is 



x,ES 



3GM sun 

^3 {Iz'z' ~ Ix'x') COS 9 



(15) 



For the non-magnetic CV DN system, M = Mi and is not in solar masses. For the misaligned CV 
DN system, the component of average torque along the line of nodes is 



x,DN 



2rf3 {Iz'z' ~ Ix'x') COS 9 



(16) 



Note that the direction is negative as indicated by z' x z. 
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2.5. Retrograde Precession 

A rotating body responds to a torque that has been applied over a small period of time by 
gaining a small amount of angular momentum, T = ^L. The direction of the small gain in 
angular momentum is parallel to the the applied torque and perpendicular to the spin axis z' of the 
body. The rotating body responds over this same small time interval by moving its spin angular 
momentum vector to a new position. By repeating this process many times over a significantly 
longer period time, the spin axis of the rotating body z' traces out a surface of a cone whose axis 
coincides with the normal to z. The direction of the trace is retrograde to the direction of the spin 
of the body and orbit as shown in Figure 1. The angular speed of the trace is the magnitude of 
the precessional angular velocity. Thus, to find the retrograde precession we need to find the time 
rate of change of the angular momentum and equate that to the torque found in Equation (14). 

To establish the total angular momentum, we need to find the total angular velocity. By 
definition, the total angular velocity of a body rotating about its axis fl is the sum of the spin rate 
N, which is along the z' axis, and the motion of its spin axis z', which is along the x' axis, 



n = Nz' + z'x— . (17) 




In this equation, we assume that no net torque is applied and l x < x > = I y i y > . From Euler's equations 
of motion, we find 

N = Iz ' z ' ~ Ix ' x ' e. (18) 
Upon substitution into Equation (1), the total angular momentum becomes 

/ „ dz'\ 

L = I z , z ,N2l + l x , x , z'x— , (19) 



dt 



of which the time rate of change is 



d r dz' T ~ d 2 z'\ . . 

dt = I ^ N ^ +I *'*\^)- (20) 

In obtaining this last equation, we assume a constant spin rate and constant moment of inertia 
coefficients. 

The time rate of change of the angular momentum is just the torque we found in the previous 
subsection. Even though the magnitude of this torque is not constant over one orbit and therefore 
the motion of the angular momentum vector moving to its new position pulses with this driving 
torque, we can assume that over long periods of time this net retrograde precessional angular 
velocity 7 is steady. Other than what is mentioned here, we do not consider nodding motions any 
further in this work. 
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After setting the time rate of change of the angular momentum equal to the torque found in 
the previous subsection and after simplifying, we obtain the quadratic 



7 2 



h'z'N 



whose roots are real and can be summed, the smaller root being 



3GM (J- v - /_,_,) , N 

+ ~2cp~ h = ° <21> 



'■'>('• M (Iz'z' -Ix'x' \ ,, 



Like for Equations (15) and (16), we can substitute M = M sun in Equation (22) for the Sun that is 
located at a distance d from the origin and M = M moon for the Moon that is located at a distance 
a from the origin 

3G ( M sun M moon \ / I z / Z i — I x i x i \ 

i = " aw (— + — ) ( w ) cos (23> 

This equation is the retrograde precession of the zero-point of Right Ascension (a.k.a., the First 
Point of Aries) in the Earth-Moon-Sun system (see any fundamentals of geophysics or celestial 
mechanics book such as Danby 1962 and Stacey 1977). 

Likewise, we could substitute M = M2 for a secondary that is located at a distance d from the 
origin 

3GM 2 ( I z ' z > - I x i x > 



This equation is the retrograde precession of a point on the rim of the spinning, mis-aligned disk. 
This point can be considered equivalent to the zero-point of Right Ascension (a.k.a., the First Point 
of Aries) in the Earth-Moon-Sun system. 

Because the disk is spinning and because it is tilted, forces that arise from the gradient of 
the potential are not colinear, and these non-colinear forces cause a net torque on the disk, one 
component of which is along the line of nodes. This component of net torque causes the spin axis 
of the disk to precess in a direction that is retrograde to the orbital direction of the secondary. The 
net motion traces out a cone whose axis coincides with the normal to the orbital plane. We note 
that this retrograde precessional equation for a spinning tilted accretion disk in a non-magnetic 
system is generic in that we have yet to take into account the moment of inertia coefficients, a 
differentially rotating disk, and any effects due to the gas stream on the accretion disk. 



3. Retrograde Precession in Tilted Disks - Solutions To Theoretical Expressions 

Retrograde precession expressions that have been found by others for various systems can also 
be found by the technique developed in §2, and we show this in this section. Insight is gained by 
comparing derivation results. In all examples, we assume Equation (24) is valid. In addition, we 
assume the disk is of constant density, orbits are circular, and Equation (3) is valid for finding the 
moment of inertia. 
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3.1. Rigid, Continuous, Circular, Tilted Disk 

We solve Equation (3) in cylindrical coordinates to find the moment of inertia of a rigid, 
continuous disk that has constant density. We also assume the radius of the disk is much larger 
than the height of the disk, r d » h c . For a thin circular disk, 



\M d r\ 
\M d r\ 



\M d r\ 



(25) 



We assign I x i x i, I y > y >, and I z > z i to the diagonal elements in the matrix, noting that I x 'x'=Iy'y' and 
thus N = e. Substituting these into Equations (24), the retrograde precession becomes 

7 = "4T^ cos (26) 

Equation (26) is the same result as that obtained by Katz et al. (1982) who approximated a disk 
as a thin ring to find the nodding motion of accreting rings and disks in X-ray binaries. In their 
work, they consider the ring/disk to be tilted and acted upon by a secondary mass. The Keplerian 
orbital angular velocity, 

lGM 2 

"-V-sr- (27> 

is easily identified in Equation (26). We assume the angular velocity of the disk around the z' axis 
is uj d = e. Thus, the retrograde precession of a thin, circular, continuous, constant density, rigid, 
tilted disk is 

3 to 2 

7 = --— cos6>. (28) 

At the edge of the continuous rigid disk, the centripetal force equals the gravitational force. There- 
fore, we can also represent the angular velocity of the disk at the disk edge by 



GMi 

^ = y— p ( 29 ) 

We can substitute this last equation and Equation (27) into Equation (28) to obtain 

7 = -^(^) 3 ^cos0 (30) 

where mass ratio q = M2/M1. 

If we assume n = 2 in the polytropic equation of state used by Larwood (1997), then their 
expression for retrograde precession for protoplanetary accretion disks reduces to our Equation (30). 
Likewise, if we assume n = 1.5 in the polytropic equation of state used by Romero et al. (2000), 
then the expression they assume for retrograde precession in quasars and black hole accretion disks 
reduces to our equation. This equation is proportional to the result obtained by Kondo, Wolff, & 
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van Flandern (1983), and subsequently by Barrett, O'Donoghue, & Warner (1988), who establish 
the precession rate of an orbiting particle in a tilted disk. Their proportionality constant is 1/2, 
however. 

Substituting Newton's version of Kepler's Third Law for the binary orbit into Equation (30), 
we obtain yet another version of retrograde precession 

7 = — - cos CM — I LO. (31) 



If we assume a small angle approximation, then Equation (31) reduces to a later derivation by 
Warner (1995, 2003) that is based on a retrograde precession of a tilted disk's outer annulus 
(Kumar 1986, 1989). The absolute magnitude of Equation (31) is the same result as that obtained 
by Osaki (1985) for prograde precession in non-tilted or non-warped disks, a result that is too high 
for prograde precession as pressure effects within the disk have been neglected (Lubow 1992). 

The expressions for retrograde precession developed by others, shown above, are the equivalent 
to our derivations if we assume that the disk rigidly rotates at the same rate as the primary, i.e., 
Cod = e. In other words, differential rotation has been neglected. 



3.2. Differentially Rotating, Circular, Tilted Disk 

In obtaining / in Equation (25), we assume the disk to be circular, thin, of constant density, 
rigid, and continuous. However, Papaloizou & Terquem (1995) note that accretion disks are not 
solid objects, and thus angular frequency is not constant with radius in fluid disks. That is, disks are 
fluid and thus individual rings can rotate differentially and this should be considered in evaluating 
Equation (25). 

From Equation (25), we note that I z i z i — I x / X r = I x i x i. We assume this same equation holds in 
a differentially rotating disk. Equation (24) simplifies to 

3GI X , X ,M 2 

7 = —- — - ^-coscV. (32) 

' 2 e I z , z , d 3 v ' 

Angular frequency is not expected to be constant with disk radius and thus Equation (29) is revised 
to its most general form, 

where r is the radial distance from the origin to any particle within the disk. As we still assume 
the disk is circular, we still expect I x 'x ,= Iy>y' even though the disk is differentially rotating. We do 
expect that the differential angular velocity to affect the I z t z i coefficient. With these assumptions 
and allowing e = u) r , the retrograde precession of a thin, circular, constant density, differentially 
rotating, tilted disk is 

3GM 2 Wcosfl 
7 " 2 d3 G i/2 M i/2 Jz>z/ V > 
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where 

rh c /2 p2-k rr d 

J 2V = / / / pr- 3/2 r 3 drd<j)dz. (35) 



-h c /2J0 JO 

Note that because the units of I z i z i change when we include the r 3 / 2 term inside the integral, 
the integral is redefined to J z ' z >. The solution to this integral is, assuming a constant den- 

5 12 

sity disk, J z t z > = (A/h)irh c rJ p. Similar calculations of J x i x > and J y i y i show that our assump- 
tion of I x i x i = Iy/yi is correct in Equation (32). As p = MdH^r^h^ 1 for a cylinder, then 
Jz'z' = (4/5) M^r^ 2 . Substituting this result into Equation (34) and substituting in our pre- 
vously found solution for I x i x i = \j\M^r\^ the retrograde precession of a thin, circular, constant 
density, differentially rotating, tilted accretion disk is 

15 M 2 / 7-^3/2 f G M 1 \ 1/2 



32 Mi V d J \ d 3 
15^2 /^\ 3 
32 Mi \ d) 



) d; d cos9 (37) 



15 q n /m3/2. 



32 v / + 9) 



cos 



To obtain Equation (37), we have substituted the angular velocity at the disk edge, Equation (29), 
into Equation (36). To obtain Equation (38), we have substituted in Newton's version of Kepler's 
third law. As we noted earlier, this derivation assumes constant surface density and constant disk 
radius r&. 

Equation (37) is the retrograde precession found by Larwood et al. (1996) for pre-main se- 
quence binary systems where they assume a polytropic gamma of 5/3 (i.e., isentropic) in their 
equation of state. If viscous dissipation occurs, then an efficient cooling mechanism is assumed 
if fluids are to remain isentropic throughout. Equation (37) is the same result as that found by 
Terquem, Papaloizou, & Nelson (1999) for protostellar disks that have tidally induced warps. Equa- 
tion (37) is proportional to that found by Larwood (1998) for forced precession in X-ray binaries. 
Larwood (1998) found a 3/7 proportionality to our 15/32 as they averaged the tidal torque over 
the disk. Equation (38) is from Montgomery (2004). 

We note that Equation (35) is not entirely correct. In this equation, we integrate from the 
origin to the outer radius of the disk. The inner rim of the disk does not extend to the origin. 
However, the outer radius of the disk is so large compared to the radial extent of the white dwarf 
or to the radial extent to the inner rim of the disk that the effect on precession is minimal. Hence, 
Equation (35) is valid for our study of precession in CV DN systems. 

In the derivation of Equations (36)-(38), we assume that the disk differentially rotates. Implied 
in this derivation is that the disk precesses as a unit. Support for this assumption is given by 
Papaloizou & Terquem (1995) who analytically show, and Larwood et al. (1996) numerically 
verify, that rigid body precession of a Keplerian disk is possible so long as the sound crossing time 
scale in the disk is small compared to the precession time scale of the disk. Further, numerical 
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simulations by Murray (1998) and observations by Patterson et al. (1998) do not indicate that 
disks differentially precess. 



3.3. Differentially Rotating, Elliptical, Tilted Disk 

In §2.4, we assume that the major components in accreting CV DN systems are the primary, 
secondary, and the accretion disk. We assume that the orbit is circular, symmetry is about the 
z' rotation axis, and the disk is circular. A circular disk is mirror symmetric and rotationally 
symmetric around its z' axis and, as a result, we could find equal body moment of inertia coefficients 
Ix'x' = ly'y' an d a constant rotation around the normal axis to the plane of the disk, Equation (24). 

However, accretion disks in CV DN systems have long been thought to be non-circular (see 
e.g., Warner 1995, 2003). An elliptical disk does not have equal radii. Thus, the moment of inertia 
coefficients about the principle axes are not equal, I x i x i / I y > y >, and T z / 0. Yet, over a time period 
of one orbit, f x still dominates. Equations (31) and (38) are still valid in non-circular disks so long 
as the disk precesses as a unit, the distance d is much greater than the radius of the primary as 
discussed in §2.2, and precession time is much longer than the orbital time. 



3.4. Rotating, Circular, Tilted Ring 

In the derivation of Equations (36)- (38), we assume that the disk differentially rotates yet the 
disk precesses as a unit. We have justified the geometry of the disk in neglecting a small, centrally 
located inner hole where the primary is located and in neglecting the outer shape of the disk being 
more elliptical than circular in CV DN non-magnetic systems. We have taken into account major 
components of the system - the primary, secondary, and the accretion disk. However, we have 
neglected the effect of the gas stream on the accretion disk. Therefore, we introduce additional 
geometrical shapes that may be used to help describe the overall geometrical shape of the disk as 
the overall geometrical disk shape may have changed due to the gas stream striking the disk. 

In this subsection, we consider the geometrical shape of a ring that is centered around the 
primary and that is tilted at the same angle as the disk. The inner rim radius of the ring r inner is 
located near the innermost rim of the disk and the outer rim radius of the ring is located at r outer - 
In this derivation, we also assume rotation and precession as a unit. 

In modeling black holes, Bardeen & Petterson (1975) state that the combined effects due to 
precession and viscosity should be to align the rotation of the inner disk annuli with the spin axis 
of the hole. Kumar and Pringle (1985) show that the alignment of the disk and the hole does occur 
in slightly tilted, steady state disks. Borrowing these ideas, we return to Equation (32) and assume 

that 6 — dinner 

3 G I X , X ,M 2 fl ,_ Q . 
7 = ---. - — -«-cos0. (39) 

^ dinner ^z'z 1 " 
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where 



GM 1 

dinner — \ ~ (40) 

V r inner 



is the angular velocity at ri nner . Note that we assume that the inner edge of the ring rotates at the 
same rate as the primary. Using Kepler's Third Law, we can relate the radius and angular velocity 
of the inner ring edge to that of the outer ring edge 

/ uj t \ 2 ^ 3 

Tinner = T uter I ~ ) (^1) 

\ dinner J 

to obtain 

M ( Tinner \ 3 / 2 h'x> (GM,) 1 / 2 

Substituting in Kepler's Third Law for the orbital period of the binary orbit, Equation (41), and 
the moment of inertia coefficients from Equation (25), Equation (42) becomes 

3 q (TouterX^I 2 . I 'outer \ , AO , 

7 = — - cos — - — uj (43) 

4 y/1 + q \ d J \Ui nner J 



•J / Tinner \ n / a a\ 

J = --q [—J— ) dinner COS 9 (44) 

3 q (Tinner^/ 2 . ..^ 

7 = --^==cos0 — — u. (45) 



or equivalently 
or equivalently 

4 ^/TTq ™" " V d 

Notice the similarities between Equations (30) and (44) and between Equations (31) and (45). 

To reduce Equation (45) further, we need to know the radius to the inner rim of the circular 
ring. The inner rim radius of the ring is unlikely to be the radius of the white dwarf as the inner 
rim of the ring is unlikely to extend down to the surface of the primary. Also, the inner rim radius 
of the ring is not likely to be the distance of closest approach of a free-falling gas stream to the 
primary (Lubow k, Shu, 1975) 

r min = 0.0488dg- a464 , (46) 

which is good for (0.03 < q < 1). Although this distance of closest approach yields a radius that is 
larger than the radius of the white dwarf, this radius is not large enough to explain a circularized 
ring: The initial path of the free-falling gas stream forms a parabolic shape around the primary, 
and after interaction with the primary's Roche lobe, the path then forms a rosette pattern. Upon 
collisions of the gas stream with itself, the rosette pattern of the gas stream settles into a circularized 
ring that has a radius that is larger than the distance of closest approach of the initial free-falling 
gas stream. Yet the inner rim radius of the ring cannot be the circularization radius (Hessman & 
Hopp, 1990) 

r circ = 0.0859d^ a426 , (47) 
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which is good for (0.05 < q < 1), as this equation describes the minimum radius of the ring's outer 
radius. 

As the ring is thin, we shall assume that r inner is half way between the minimum radius and 
the circularization radius 

/ r min ~\~ r circ \ r ao\ 

rinner = I ^ I • (48) 

As Warner (2003) finds that 

3 q ( 1.38r min \ 3/2 . 
7 = 7 n — = cos " ; w 

= _^_^ cos ^( .0673r°- 464 ) 3/2 ^. (49) 

In summary, we assume a thin ring that is centered about the origin and is located near the inner 
edge of the disk. This ring retrogradely precesses as a unit. We note that we could have repeated 
the process of differential rotation that we introduced in §3.2. However, because the ring is thin, 
we find that Equation (49) is minimally affected. 



3.5. Negative Superhump Period Excesses 



As the orbital period is slightly longer than the period of the secondary encountering the same 
disk node per orbit, we can relate the negative superhump period P and the orbital period to the 
retrograde precession period by 



r r --r- ~^orb- 

Another relation is the negative superhump period excess 



e_ = 1- 



Porb 



1 + 



P, 



orb 



(50) 

(51) 
(52) 



where we substitute Equation (50) into Equation (51) to obtain Equation (52). For example, we 
can find the negative superhump period excess for the rigid, continuous, circular, tilted disk from 
§3.1. If we assume that the radius of the disk is that given by Paczynski (1977) 

0.6d 



1 + q' 



(53) 



which is valid for mass ratios in the range 0.03 < q < 1, then we can combine Equations (31) and 
(52) to obtain 

-l 



rigid 



1 



1 + 



0.349g cos i 
(l + q) 2 



(54) 
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Likewise, we can find the negative superhump period excess for the differentially rotating, 
circular, tilted disk from §3.2 by substituting Equation (53) into Equation (38) and then into 
Equation (52) to obtain 



€—Aisk 



1 



1 + 



0.218<? cos6> 
(1 + q) 2 . 



(55) 



Likewise, we can find the negative superhump period excess for the rotating, circular, tilted 
ring from §3.3 by substituting Equation (49) into Equation (53) and then into Equation (52) to 
obtain 

-l 

(56) 

Note that we have dropped the negative sign that redundantly indicates retrograde precession. 

Equation (55) appears in Montgomery (2004). If we assume a primary mass Mi=0.8M & and 
mass ratio q = 0.4, then the secondary mass M2=O.32M . The radius of the disk is ~ 0.43<i, 
where d ~ 1.23i?o (Montgomery 2009). The rotation period at the disk edge, determined from 
Equation (29), is around one hundred seconds. This rate is fast compared with the orbital period 
Porb = 3.57 hours found using Newton's version of Kepler's third law. From Equation (55), we find 
P r ~ 22.5P or b =80 hours, assuming small angle approximation for disk tilt. 

Instead of mass ratio, we can express Equations (54), (55), and (56) in terms of the orbital 
period. For this substitution, we shall assume non-magnetic, accreting binaries that follow the 
Smith &; Dhillon (1998) secondary mass-period relation, 

M 2 /M Q = (0.038 ± O.OOS)^ 8 ±a09 (57) 

where P or b is in hours. We do not show these lengthy negative superhump period excess equations 
in this work, but they can easily be found by variable substitution where q = ^ and Mi = O.8M . 



,rmg 



1 + 



Igcos^ 



0.464 ~\ 3 / 2 



(0.0673g- a464 ) 



4. Numerically Simulated Retrograde Precessions 

In Montgomery (2009), we simulate non- magnetic CV DN systems that have mass ratios within 
the range (0.35 < q < 0.55) and low mass transfer rates. By limiting our mass ratio range to longer 
orbital periods, we eliminate systems that apsidally precess and thus show positive superhumps in 
their light curves, thereby minimizing variables. We adopt a Smith & Dhillon (1998) secondary 
mass-period relation which is for main-sequence dwarfs as the secondaries. Because we simulate 
accreting, non-magnetic CV DN systems in or above the period gap, we choose an average white 
dwarf primary mass M\ = O.8M0, the average value found by Smith & Dhillon (1998) for this 
mass ratio range. We refer the reader to Montgomery (2009) for details regarding the code yet we 
provide a summary here. 
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For our numerical simulations, we use Smoothed Particle Hydrodynamics (SPH) to simulate 
accretion disks that have constant and uniform mass particles. In the code, we assume assume a 
constant and uniform smoothing length h and an ideal gamma-law equation of state P = (7 — l)pu 
where P is pressure, 7 is the adiabatic index, p is density, and u is specific internal energy. The 
sound speed is c s = ^7(7 — I) - "- The code utilizes a Lattanzio et al. (1986) artificial viscosity, and 
we choose a = (3 = 0.5 viscosity coefficients. Our viscosity is approximately equivalent to a Shakura 
& Sunyaev (1973) viscosity parametrization {y = a'c s H where H is the disk scaleheight) a'=0.05. 
As Smak (1999) estimates a' ~ 0.1 - 0.2 for DN systems in high viscosity states, our simulations are 
more for quiescent systems. As only approaching particles feel the viscous force and since neither 
radiative transfer nor magnetic fields are included in this code, all energy dissipated by the artificial 
viscosity is transferred into changing the internal energies. 

As radiative cooling is not included in this code, an adiabatic index 7=1.01 is incorporated 
to prevent internal energies from becoming too large. By integrating the changes in the internal 
energies of all the particles over a specific time interval n, variations in the bolometric luminosity 
yield an approximate and artificially generated light curve 

L n = £duf. (58) 

i 

In all simulations, we build a disk to 100,000 particles. At orbit 200, we artificially rotate the 
disk out of the orbital plane 5° to induce negative superhumps in the light curve and to force 
the disk to precess in the retrograde direction. For example, if q = 0.4 then M2=0.32 Mq. Our 
simulation unit length scales to d ~ 1.23i?0 if we assume the secondary mass-radius relationship 
R 2 = (M 2 /M Q )- 13 / l5 R & (Warner 2003) that applies for 0.08 < M 2 M^ < 1.0 and the Eggleton 
(1983) volume radius of the Roche lobe secondary 

R2 = 0.49g 2 / 3 

a 0.6g 2 /3 + l n (i + g i/3)- > 

In this equation, R 2 and M 2 are in solar radii and solar mass, respectively. The Eggleton (1983) 
relation is good for all mass ratios, accurate to better than 1%. Taking the radius of the primary 
to be 6.9xl0 8 cm, the scaled radius of the white dwarf is R\ = 0.0081(1 

After building the disk, we allow the disk to evolve in the short term to a quasi-equilibrium 
state where particles are injected at the rate they are removed from the system by either being 
accreted onto the primary or the secondary or lost from the system. The average net rate of 
accretion for a quasi-static disk is around 500 particles per orbit as shown in Montgomery (2009). 
As this net rate is approximately half that injected to build the disk, the steady state mass transfer 
rate reduces to m ~ 1 x lO~ lo M0/yr or a rate similar to an SU UMa in quiescence. We estimate 
the density of the gas near L\ to be p ~ 10~ 10 g cm -3 using m p = (4/3)-7r/i 3 p. As we maintain 
100,000 particles in the steady state disk, then the total mass of the disk is M^= 100,000 x m p 
= 3.5xl0 22 g or M d ~ 2xlO" n M , a negligible value compared to the mass of either star. If we 
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assume the outer radius of the disk is Ra « 2r c j rc , or twice the circularization radius (Warner 2003) 

^ = 0.0859(?- a426 (60) 
a 

where 0.05 < q < 1, then for q = 0.4 we find r circ ~ 0.13d and R d ~ 0.25d. As p « M d R^ 3 , self 
gravity is neglected in these simulations. 

In Montgomery (2009), we find the calculated negative superhump period, its propagated error, 
and the nodal superhump period excess as output to the code. We list these values in Table 1. 
The units for period are hours. We also find the source that powers the negative superhump. Our 
results suggest that negative superhumps are due to an increase in the number of particles that 
reach innermost annuli of a tilted disk (Montgomery, 2009). As inner annuli have higher energies, 
the combination of increased number of particles and higher energy inner annuli yields more light 
being emitted. This additional light waxes and wanes with the amount of gas stream overflow and 
the amount of particle migration into innermost annuli. As the secondary orbits the tilted disk, the 
gas stream mostly oveflows one face of the tilted disk for approximately the first half of an orbit 
and then over the other face of the tilted disk for approximately the second half of an orbit. As a 
result of the disk tilt and thus gas stream overflow and additional particle migration to innermost 
annuli, a yellow colored ring that is denser that its surroundings is generated near the innermost 
annuli (see Figure 3). This denser ring is important to the generation of the negative superhump. 

We do note that the relative surface densities across the face of the disk vary by only ~ 4% 
from annuli to annuli, regardless if the disk is tilted or not. Exceptions are the disk rim and a denser 
ring, colored green, located near ~ 0.2d in Figure 3. The latter is caused by gas stream overflowing 
the tilted disk and striking the disk face. Of the two accretion disks shown in Figure 3, the tilted 
disk is more blended in density. Hence, assuming a constant density disk in our derivations of 
retrograde precession in §3 is justified. 

Note in Figure 3 that the ring is tilted at the same angle as the disk and that the ring is 
rotating. Therefore, in addition to the disk, the ring is also subject to precession. Also note that 
numerically simulated CV DN disks (Figure 3) are elliptical and do not have excessively large, 
centrally located inner holes. As the disk radius is much larger than the radius of this small hole, 
our previous assumption of integration from the origin of the system to the outer edge of the disk 
is justified. 



5. OBSERVATIONAL DATA 

Table 2 lists observational data including orbital periods, negative superhump periods, and 
negative superhump period excesses for several non-magnetic, accreting CV DN systems. If errors 
are known, they are listed in parentheses. In the table, NL is novalike, PS is the permanent 
superhumper category of NL, IP is intermediate polar, SW Sex and VY Scl are members of NLs, 
ER UMa is a member of the SU UMa subclass. 
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In Table 2, we also list positive superhump periods and precessional data for those systems 
that precess both in the prograde and in the retrograde direction. Prograde precession is due to 
density waves forming radially through the middle-to-outer annuli of the disk (Smith et al. 2007) 
and is described by a dynamical prograde term plus a retrograde pressure term (Lubow 1991a,b). 
The retrograde precession term is a correction to the prograde precession term due to these spiral 
density waves. In this work, we do not establish the geometrical shape of an accretion disk with 
these radial density waves as these radial density waves are not generated in long orbital period 
CV DN systems. We leave this as future work. 



Figure 4 shows a comparison of the observational nodal superhump period excess (□ for IPs, 
o for AM CVn, * for CV DN that progradely and retrogradely precess, and X for those systems 
that only retrogradely precess) to the numerical negative superhump period excess (A) values 
listed in Tables 1 and 2, respectively, as a function of orbital period. The width of the symbol 
can be taken as the error in the orbital period. In addition, we overplot a short dashed line that 
represents a retrogradely precessing rigid disk - Equation (54), a dash-dot line that represents 
a retrogradely precessing but differentially rotating disk - Equation (55), a long dash line that 
represents a retrogradely precessing ring - Equation (56), and a solid line that represents a disk 
that is differentially rotating yet retrogradely precessing as a unit and has a precessing ring. 

As shown in Figure 4, the retrograde precession of an inner annuli ring (long dashed line) does 
not match any of the observations, a result in agreement with the idea that the disk cannot be 
represented by a ring. The retrograde precession of solid, rigid disk also does not match any of the 
observations, a result in agreement with the idea that the disk rotates differentially as the disk is 
a fluid. 

Also shown in Figure 4, the retrograde precession of a differentially rotating disk does not 
match the numerical simulations. Of the four lines shown in Figure 4, the numerical simulations 
are best described by the solid line. To generate the solid line in Figure 4, we assume that the net 
retrograde precession is due to a differentially rotating disk that has superimposed on it a rotating 
ring due to the gas stream overflowing the tilted disk and particle migration as the secondary orbits. 
The net precession is 



as shown in Figure 4. 

Smith & Dhillon (1998) find that CV secondary stars with orbital periods shorter than 7-8 
hours are indistinguishable from main-squence stars in detached binaries in terms of spectral type, 



6. 
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mass, and radius. They find that CVs with periods up to 7-8 hours almost always have main- 
sequence secondaries, although they do point out that systems with peculiar secondaries do exist. 
We assume their secondary mass-period relation in both our theoretical expressions and in our 
numerical simulations. If Smith & Dhillon (1998) and this work are correct, then the solid line in 
Figure 4 should represent those systems that have main sequence secondaries. The one observational 
system that is on the solid line is VI 193 Ori. HS 1813+6122 is nearby to one of our numerical 
simulation points. As Smith & Dhillon (1998) find that the average white dwarf mass above the 
period gap is M\ = 0.80 ± O.22M0, the same average mass we chose for the numerical simulations 
and theoretical expression for retrograde precession, we find M 2 ~ 0.267M Q and q ~ 0.33 for V1153 
Ori. We find M 2 ~ 0.281M© and q ~ 0.35 for HS 1813+6122. In agreement with the ideas of Smith 
& Dhillon (1998), these above-the-period-gap systems have stellar masses. 

We do note the large scatter in the observational data. As shown in Figure 4 by the various 
lines, the geometric shape of the disk significantly affects the precession rate and could be a major 
source to the scatter. Another smaller source to the scatter not considered in this work is density 
waves that may be generated in the disk. Therefore, we have identified the observational data in 
Figure 4 by various different symbols to indicate systems with different geometries or otherwise 
unusual disks. For example, those systems in Table 2 labeled IPs are expected to have a large 
absence of inner annuli due to the primary's large magnetic field that disrupts the disk. In this 
work, we did not find the retrograde precession of a differentially rotating, retrogradely precessing, 
hollow cylinder and thus we flagged these observational systems in Figure 4 as squares. A hole as 
large as that experienced by an IP's disk would not be described by the solid line or any other 
line in Figure 4 as the geometrical shape of this system is different than that described by these 
lines. Likewise, the geometrical shape of an accretion disk that both progradely and retrogradely 
precesses would not be described by any of the lines in Figure 4. These systems have radial density 
waves in their accretion disks and may also have an innermost ring due to gas stream overflow 
and particle migration to innermost annuli because the disk is tilted. We did not consider these 
geometrical shapes or disks with density waves in this work. We leave these different geometrical 
shape disks for future work. 



7. DISCUSSION 
7.1. Warp or Tilt? 

Warps in disks can be localized, say in inner regions or in outer regions (e.g., Boffin et al. 
2003), or across the entire disk. An example of localized warps are those generated in inner 
disk annuli. Nelson & Papaloizou (1999) find inner disk warps that are induced by misaligned 
angular momentum vectors of the disk and the rotating black hole. Petterson (1977) and Iping 
and Petterson (1990) study radiation pressure as a source to disk warp. If a primary is a source 
of strong radiation, then the surface of the accretion disk can absorb this radiation and re-emit it 
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normal to the disk surface, inducing a torque on the disk that may warp the disk if the central 
source luminosity exceeds a critical value (see e.g., Petterson 1977; Pringle 1996, 1997; Wijers & 
Pringle 1999). The strong radiation source twists the disk so that each annuli affected precesses at 
the same rate. In these studies, the warp precesses relative to the inertial frame as a rigid body. 
Bisikalo et al. (2004) study CV accretion disks and they too find a localized inner disk warp that 
precesses yet the disk as a whole does not precess. 

An example of entire disk warp that precesses as a unit are the studies of Larwood et al. (1996), 
Larwood & Papaloizou (1997), and Larwood (1997). These authors find precession rates due to 
a tidal interaction between a circumprimary disk or a circumbinary disk and a secondary point 
mass on a parabolic trajectory in protostellar systems. Their results indicate that disks become 
modestly warped (1° — 6°) during non-coplaner encounters and that disks precess approximately as 
a unit. Murray et al. (2002) find that a warp can be induced in CV accretion disks by an inclined 
magnetic dipole field that is centered on the secondary. They also find that their warp precesses 
retrogradely and as a unit. 

An example of a localized warp that may cause the entire disk to warp is the studies by Foukes, 
Haswell, & Murray (2006) who generate warps in X-ray binary disks. Their source to warp is also 
radiation driven by the primary. They find that inner disk annuli can warp and the warp precesses 
nearly as a rigid body in the retrograde direction. For less extreme mass ratio systems, they find 
that inner and outer regions of the disk warp to the point that they can say that the entire disk 
is tilted to the binary plane. As a result of the entire disk tilt, the whole disk precesses in the 
retrograde direction, a result also found in Montgomery (2004, 2009). 

For CV DN systems, warps seem unlikely and are usually not considered (e.g., Boffin et al. 
2003) as CV DN primary are not known to be strong radiation sources. In addition, this study 
considers only CV DN that are considered non-magnetic. The geometry suggested in Larwood et 
al. (1996), Larwood & Papaloizou (1997), and Larwood (1997) is unlikely for CV DN systems. 
Also unlikely is a misaligned disk and rotating primary angular momentum vectors as this source 
seems to only warp inner annuli without the generation of negative superhumps (Foulkes, Haswell, 
& Murray 2006). Disk tilts have been postulated to involve the entire disk or only a portion such 
as the outer edge (Bonnet-Bidaud, Motch, & Mouchet 1985; Kumar 1986, 1989). As Foulkes, 
Haswell, & Murray (2006) and Montgomery (2009) find that whole disk needs to be tilted in order 
for the gas stream to flow over and under the disk edge, entire disk tilts are the likely geometry 
for non-magnetic, accreting, tilted, spinning systems if negative superhumps of the correct shape 
and frequency per orbit are expected to be seen in the light curves. The source to disk tilt in 
non-magnetic systems remains unknown, and we leave this as future work. As Figure 5 shows, a 
fully tilted disk may look like a disk warp as the system retrogradely precesses. In Figure 5, we 
only show one half of the total disk precessional period. 
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7.2. Whole or Partial Disk Precession? 

In Montgomery (2009), we find that if the entire disk is tilted, then the gas stream can overflow 
a disk face and reach inner annuli. As the secondary orbits, the gas stream flows over one face of the 
disk for ~ 1/2 an orbit and the flows over the other face of the disk for ~l/2 an orbit. The location 
where the gas stream transitions on the disk rim from flowing over to under and vice versa (i.e., the 
nodes) precesses in the retrograde direction. Like the slow Westward movement of the First Point 
of Aries along the Equator of the Earth, the line of nodes also has a slow Westward movement 
along the rim of the disk. If the disk precesses as a unit and the line of nodes shifts Westward by 
a small amount each orbit, then the location where the gas stream strikes the innermost annuli 
must also shift Westward by a small amount each orbit. That is, the spinning, tilted ring formed 
near the innermost annuli of the disk and formed by the gas stream overflowing the tilted disk also 
precesses in the retrograde direction. Therefore, to explain precession in our numerical simulations, 
we net the retrograde precession of the tilted, rotating ring and the tilted, rotating disk. 

Bisikalo et al. (2004) finds that only their spiral density wave retrogradely precesses, not 
the entire numerical simulated CV accretion disk. Observations by Patterson et al. (1998) do 
not indicate that disks precess differentially Numerical simulations by Murray (1998) do not 
indicate that disks precess differentially. Our results do not indicate that disks precess differentially. 
Papaloizou & Terquem (1995) analytically show, and Larwood et al. (1996) numerically verify, that 
rigid body precession of a Keplerian disk is possible so long as the sound crossing time scale in the 
disk is small compared to the precession time scale of the disk. We find in Montgomery (2009) 
that the entire disk needs to tilt for additional particle migration to reach inner annuli and produce 
negative superhumps. We find in this work that both outer and inner annuli precess in fully tilted 
disks. These results suggest that partially retrogradely precessing non-magnetic CV DN tilted 
accretion disks are unlikely. 

7.3. Precession Source 

Driven precession involves torques on a body. One example is a wobbling coin on a table. The 
coin spin rate is slow but its wobble rate is fast due to a torque on the body. Another example is a 
spinning top inclined at an angle from the vertical. The spin rate is fast but the precession rate is 
slow due to torques on the body. The precession of the Earth due to tidal torques by the Moon and 
the Sun is another example of a slow precession rate compared to Earth's (relatively) fast spin rate. 
For a CV DN system with q = 0.4, we previously calculated the retrograde precessional period 
to be approximately 80 hours. This period is slow compared to the primary's rotation period of 
around one hundred seconds, if we assume that the primary's rotation period is similar to that at 
the disk edge. Therefore, driven precession could be a viable source of retrograde precession in CV 
DN systems that show negative superhumps in their light curves. Katz (1973), Roberts (1974), 
Katz et al. (1982), Kondo et al. (1983), Kumar (1986, 1989), Barrett et al. (1988), Warner (1995, 
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2003) are all examples of those employing or assuming torques to induce retrograde precession. 

As we show in Figure 4, retrograde precession in CV DN disks is likely due to the same 
source that causes the Earth to retrogradely precess, that is, by tidal torques. The difference in 
precessional values of the Earths' equinoxes and a CV DN's line of nodes is due to the geometric 
shape of the precessing object (e.g., oblate spheroid, disk, etc.), the state of matter (i.e., fluid, solid, 
etc.), and any effects due to the accretion stream. 

Bisikalo et al. (2004) generate a spiral density wave in innermost annuli of their numerical 
simulated CV disks. Their spiral density wave is caused by sheared elliptical inner disk annuli. 
They find that only the spiral density wave retrogradely precesses, and they cite Kumar (1986) as 
the source to retrograde precession. However Kumar (1986) finds the retrograde precession due 
to a tilted disk's outer annuli and thus cannot be the source to retrograde precession of the spiral 
density wave. However, we may be able to help Bisikalo et al. (2004). Lubow & Shu (1975) find 
that pressure can be generated in the disk due to shocks that are induced by gas stream overflow, 
and Lubow (1991a,b) finds that additional pressure due to spiral density waves could cause the 
disk to retrogradely precess. If the inner disk has a spiral density wave instead of a dense ring, then 
the net retrograde precession can still be explained by the solid line in Figure 4. That is, the disk 
would retrogradely precess due to combined effects due to tidal torques and pressure from a spiral 
density wave located near the innermost disk annuli (instead of a ring). We plan to show this in a 
future work. 



8. SUMMARY, CONCLUSIONS, AND FUTURE WORK 

Negative superhumps are found in observational systems with high mass ratios, long orbital 
periods, and high mass transfer rates. In Montgomery (2009), we find that the negative superhump 
can be due to extra light being emitted from inner disk annuli. The extra light is modulated by gas 
stream overflow that strikes inner disk annuli and by particle migration as the secondary orbits. 
For this result to be true, we find in Montgomery (2009) that the disk must be fully tilted. The 
source to disk tilt has not yet been resolved and we leave this as future work. Fully tilted disks 
could be subject to driven torques that cause the disk to retrogradely precess. 

Patterson et al. (1993) suggest that the retrograde precesison in CV DN accretion disks may 
be by the same source that causes the Earth to retrogradely precess. Using the Earth-Moon-Sun 
system as a model, we generate a generic theoretical expression for a spinning disk that is misaligned 
with the orbital plane to test this idea. We assume the potential of MacCullogh is valid for our 
spining, tilted disk that contains a massive primary that is smaller than the thickness of the disk 
and the distance to the secondary. Our expression matches that which describes the precession of 
Earths' equinoxes. 

By making appropriate assumptions, we reduce our generic expression for retrograde preces- 
sion to those generated by others, or to those used by others, to describe retrograde precession in 
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protostellar, protoplanetary, X-ray binary, non-magnetic accreting CV DN, quasar and black hole 
systems. Because our generic expression reduces to those generated by others for a variety of sys- 
tems, retrograde precession by tidal torques on a spinning, misaligned primary may be ubiquitous. 
Previous expressions treated the disk or ring as rigid and/or didn't include the effects on the disk 
by the accretion stream. Because the disk is fully tilted, the gas stream overflows one face of the 
disk and particles migrate to inner annuli, creating a spinning, tilted ring that is denser than its 
surroundings as shown in our numerical simulations. 

In this work, we generate a theoretical expression for the retrograde precession of a spinning, 
tilted ring. We assume that the rotation rate of the primary matches that of the inner rim of the 
spinning, tilted ring. This implies that the axes of the ring and the primary align. We show that 
the net retrograde precession in CV DN non-magnetic accretion disks needs to include the effects 
due to both the precessing disk and the precessing ring. If we include the net moment of inertia 
due to the disk and the ring, then our theoretical expression agrees well with our numerical results. 
Our findings support the suggestion by Patterson et al. (1993) that tidal torques like those by the 
Moon and the Sun on the spinning, tilted Earth are the source to retrograde precession in spinning, 
tilted CV DN accretion disks. 

In our theoretical expression and our numerical simulations, we assume an average primary 
mass Mi = O.8OM that matches the average primary mass Mi = 0.80 ± 0.22M Q found by Smith 
& Dhillon (1998). If the results in Smith & Dhillon (1998) and this work are correct, then from 
Figure 4 we find that V1193 Ori and HS 1813+6122 secondaries may have stellar masses. 

The scatter in observational precessional values seems to be mostly from differing geometrical 
disk shapes. The scatter could also be from any density waves that may be present in the disk. 
Bisikalo et al. (2004) finds that a spiral density wave is generated in innermost annuli of their 
numerically simulated CV accretion disks. They find that their spiral density wave retrogradely 
precesses. Pressure effects generated within the disk from a spiral density wave could add to 
retrograde precession. The disk could retrogradely precess due to combined effects due to tidal 
torques and pressure from a spiral density wave. We plan to show this in a future work, and we 
anticipate that the result is the same solid line shown in Figure 4. However, unlike that found by 
Bisikalo et al. (2004), we expect the disk to still precess as a unit. 

We do note that we did not consider significant mass loss of the secondary over time and this 
may also affect precessional values. If the mass lost by the secondary is not significantly acquired 
by the primary and is instead mostly lost via e.g., winds, jets, or circumbinary disks, then the 
geometric shape of the disk, moment of inertia, mass ratio, and/or orbital period changes and thus 
so changes the precessional values. Whether the retrograde precession and orbital period values 
change one-to-one is not fully understood and we leave this as future work as well. 
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Fig. 1. — Sketch of the Earth-Sun binary system (top panel) and of the non- magnetic CV DN 
system (bottom panel). Thin solid line represents the axes of the ecliptic with rotation in the 
prograde, counter-clockwise direction around the z-axis with speed Co. The primary (i.e., the Earth 
in the top panel or the disk it the bottom panel) is at the origin and rotates counter-clockwise with 
speed e around the dashed z' axis where (x',y',z') are the body axes. The z' axis is at an obliquity 
angle 6 relative to the z-axis. Note that the primary in the CV DN system (bottom panel) is small 
compared to the size of the disk and is thus hidden from view. The net torque around the z-axis 
is caused by differing gravitational attraction forces between opposite points on the surface of the 
Earth (top panel) or disk (bottom panel) and the Sun (top panel) or secondary mass (bottom 
panel), shown as a sphere on the y-axis at a distance d from the origin. The net torque causes the 
rotation axis of the primary, z', to precess in the retrograde direction. Schematics are not drawn 
to scale. 
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Fig. 2. — Schematic of a non-spherical attracting body m that has a small element dm located at 
distance r' from its center of mass and at an angle a from the line of centers connecting the center 
of mass of the attracting body with the center of mass of a far away point mass body M. The 
far-away body is located at a distance d where d» r'. Schematic is not to scale. 
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Fig. 3. — Relative densities of tilted (left panel) and non-tilted (right panel) disks . Both panels 
show orbit 220, frame (220.0) of 200 possible frames per orbit. The 20° tilted disk shows a 
green higher density ring near where the gas stream strikes the disk face. The secondary is not 
shown but the gas stream is shown leaving the inner Lagrange point. In the figure, relative den- 
sity is scaled orange/red=0.6%, yellow/red=4%, green/red=8%, white/red=18%, blue/red=29%, 
violet/red=83%. 
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Fig. 4. — Negative superhump period excess as a function of orbital period. Observational data are 
shown as stars, X's, squares, and diamonds. IPs are squares. AM CVn is a diamond. Systems that 
only show negative superhumps in their lightcurves are X's, and systems that show both negative 
and positive superhumps in their light curves are stars. Numerical data are shown as triangles. 
The short dash line represents a disk that rigidly rotates at the same rate as the primary. The solid 
line represents a disk that rotates differentially with a superimposed ring on the inner edge of the 
disk. The dash-dot line represents a disk that rotates differentially. The long dash line represents 
a rotating ring. Data are from Tables 1 and 2. 
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Fig. 5. — Retrograde precession of a disk fully tilted 5° around the line of nodes, out of the orbital 
plane, showing that a fully tilted, precessing disk can look like a warped disk. 
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Table 1: Negative Superhump Simulation Data 
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Table 2: Positive and Negative Superhump Observational Data, 
Excesses 
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0.12721(6) 


1.561 


4.306 


0.1385(2) 


0.12356 


0.0887(16) 


0.029 


DW UMa c '" 


8 


0.33 aa 
0.255 c 
> 0.24 6b 


0.13661 


2.2573 


5.038 


0.14541 


0.1330 


0.0644 (20) 


0.026 


PX And c <f ' cc 


8 


0.329(ll) c 


0.146353(1) 


1.776 


4.228 


0.1595(2) 


0.1415 


0.0898 


0.0331 


BH Lyn c ' dd 


8 


n /icee,+0.15 

u - 40 -o.io 
0.41 ± 0.26* 

0.301(15) c 


0.15575(1) 


2.128 


3.413 


0.16805 


0.1490 


0.079 


0.0433 


RX 31643" 


8 




0.120560(14) 




3.917 




0.11696(8) 




0.032(2) 


RR Cha-^ 


6 




0.1401 


4.705 


5.049 


0.1444 


0.1363 


0.031 


0.027 



Table 3: Positive and Negative Superhump Observational Data, Calculated Precessional Periods, and Period 
Excesses (Continued from Table 2) 



System 


Type 


q 


^orb{ a ) 


Ppa(d) 




Pa{d) 


Pn{d) 






AT Cnc flfl 


3,7 


0.32 - 1.04 


0.2011(6) 




5.051 
11.03 




0.1934(8) 
0.1975(8) 




0.0179(10) 
0.0383(10) 


IR Gem hh ' ii 


4 


0.154 c 


0.06840 


1.390 


2.105 


0.07194 


0.0663 


0.052 


0.031 


TX Col-" 


6 


> 0.33 


0.2375 


1.204 


1.669 


0.2958 


0.2083 


0.25 


0.1229 


SDSS J0407 fcfc 


4? 




0.17017(3) 




6.727 




0.166(1) 




0.024(1) 


SDSS J2100" 


4 




~0.0833 


1.767 


4.083 


0.0875 


0.0817 


0.050 


0.020 


CAL 86 mm 


4 




0.06613 




1.403 




0.06313 




0.045 


KR Aur" n 


2 


0.60±°° 


0.1628 




4.489 




0.1571(2) 




0.0350(2) 


HS 1813+6122PP 


8 




0.1479 




3.166 




0.1413 




0.0446 


V2574 Oph" 


9 




0.14773 




3.429 




0.14164 




0.412 


RX 1643+34" 


8 


°- 45 0.20 


0.12050(14) 




3.917 




0.11696(8) 




0.03000(8) 


AY Psc rr 


1 


0.21732(9) 








0.2057(1) 
0.2063(1) 
0.2072(1) 




0.0535 
0.0506 
0.0466 


BF Ara ** 


1 


~0.21 


0.084176(21) 






0.08797(1) 


0.082159(4) 


0.0451(3) 


0.0244(2) 


VI 193 Ori MU 


3 




0.1430 








0.1362 




0.0476 



X NL, 2 VY Scl/NL, 3 PS/NL, 4 SU UMa/DN, 5 ER UMa/DN, 6 IP/NL, 7 Z Cam, 8 SW Sex/NL, 9 Nova, 
?classification uncertain 

"Taylor et al. (1998), 6 Huber et al. (1998), c Patterson et al. (2005), d Kraicheva et al. (1999), e Wu et al. 
(2002), -^Warner (1995), ^Borczyk et al. (2003), ^Patterson et al. (1997), * Arenas (2000), ' Harvey et al. 
(1995), ^Patterson et al. (1995), 'Thorstensen et al. (1997), m Skillman et al. (1999), "Nelemans et al. 
(2001), "Roelofs et al. (2006) p Retter et al. (2003), ^Hellier & Buckley (1993), r Schoembs (1982), s Barrera 
& Vogt (1989), *Gao et al. (1999), "Patterson et al. (2001), "RX J1643. 7+3402 - Patterson et al. (2002), 
w Shafter & Szkody (1983), ^Retter et al. (1997), ^Patterson (1998), 2 Rodriguez- Gil et al. (2007b), aa Biro 
(2000), b6 Araujo-Betancor et al. (2003), cc Boffin et al. (2003), ^Dhillon et al. (1992), ee Hoard & Szkody 
(1997), ff Woudt & Warner (2002), fSKozhevnikov (2004), hh Li et al. (2004), "Fu et al. (2004), ^Retter 
et al. (2005), fefc SDSS J040714.78-064425.1 - Ak et al. (2005), "SDSS J210014. 12+004446.0 - Tramposch 
et al. (2005), mm Woudt et al. (2005), ""Kozhevnikov (2007), °°Kato et al. (2002), w Rodriguez-Gil et al. 
(2007), ^Kang et al. (2006), rr Gulsecen et al. (2009), ss Szkody & Howell (1993), "Olech et al. (2007), 
uu Xk et al. (2005b) 
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